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Electron scattering on a thin layer where the potential depends self-consistently on the wave 
function has been studied. When the amplitude of the incident wave exceeds a certain threshold, a 
soliton-shaped brightening (darkening) appears on the layer causing diffraction of the wave. Thus 
the spontaneously formed transverse pattern can be viewed as a self-induced nonlinear quantum 
screen. Attractive or repulsive nonhnearities result in different phase shifts of the wave function on 
the screen, which give rise to quite different diffraction patterns. Among others, the nonlinearity 
can cause self-focusing of the incident wave into a "beam", splitting in two "beams", single or 
double traces with suppressed reflection or transmission, etc. 
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The spontaneous formation of spatial structures (pat- 
terns) due to nonlinearity is well known far dissipa- 
tive systems driven away from equilibrium.EJ In solid 
state physics those patterns have been mostly studied i: 
the regime governed by classical macroscopic processes,! 
where quantum coherence effects were not important. In 
this paper we predict the spontaneous formation of quan- 
tum coherent non-dissipative patterns in semiconductor 
heterostructures with nonlinear properties. 

Since the Schrodinger equation is linear, the nonlinear- 
ity appears in quantum systems due to the many-body 
effects and/or the coupling with the environment. In a 
mean-field approximation this problem can be traced to 
the self-consistent Schrodinger equation with the Hamil- 



tonian H 



2m 



V 2 + V(r) + V oS [\ip(r)\ 2 ], where in ad- 



dition to the external potential V^(r) the self-consistent 
potential V e g is introduced, representing a nonlinear re- 
sponse of the medium.13 The potential V c s depends on the 
probability IVK 1 ")! 2 of the carrier to be located at r. When 
(in a weakly-nonlinear case) it is proportional to that 
probability, the resultant equation for a single-particle 
wave- function ip(r) is the so-called npnlinear Schrodinger 
equation (NSE) with a cubic termQ encountered in dif- 
ferent contexts of the solid state physics: (i) the polaron 
problem^ where the strong electron-phonon interaction 
dcformates the lattice thereby providing an.-|attractive 
potential;Q (ii) the magnetopolaron problcmEl in semi- 
magnetic semiconductors, where the exchange interac- 
tion between the carrier spin and the magnetic impuri- 
ties leads also to an effective attractive potential;Ero (iii) 
Hartree-type interaction between electrons, giving a re- 
pulsive potential,E3 and othersc 

Motivated by the great progress in heterostructure 
fabrication, some important results have been obtained 
recently in the framework of the cubic NSE for the 
situations when the nonhnearities are concentrated, in. 
thin semiconductor layers modeled by 5-potentialsJ3'M~tJ 
Among these results, we may mention the multiplicity 
of stable states found in different physical situations for 
which tunneling is important: an array of .scmimagnetic 
quantum dots,H a quantum molecular wire,EJ a doped su- 
perlattice formed by <5-barriers.ll3 Another is the oscilla- 
tory instability of the flux transmitted through the non- 



linear layer .El It should be noted, however, that all these 
results are restricted to one-dimensional spatial supports, 
which means that the longitudinal and transverse degrees 
of motion are assumed to be decoupled. Disregarding 
that assumption in this paper, we show that consider- 
ing additional spatial dimensions opens up the possibil- 
ity of qualitatively new nonlinear phenomena - sponta- 
neous formation of spatial transverse patterns, which are 
quantum- mechanically coherent. 

Consider a thin layer in the xy-plane with the con- 
centrated nonlinearity. We model the layer by using the 
(5-function which simplifies greatly the calculations with- 
out modifying qualitatively the results. Keeping in mind 
possible pattern formation and analogy with the optics, 
the layer can be thought of as a screen. The steady-state 
scattering problem for the thin ^-layer is governed by the 
NSE 



Zm 



+ {A + B\ijj{Y)\ 2 )8{z) v b{r)=Evb{v). (1) 



The external potential A is allowed to be of the both 
signs, i.e., A > if it is a barrier and A < if it is a well. 
B is the strength of the nonlinear potential: B < for the 
attractive and B > for the repulsive interaction. We 
do not specify the concrete physical model, because our 
results could be applicable to any of the above mentioned 
systems, although the most feasible candidates for the at- 
tractive case are believed to be semimagnetic heterostruc- 
tures like CdTe/CdMnTe and CdTe/HgCdMnTe, where 
both, the height A of the barrier and the strength B can 
be varied by choosing the alloy composition. In addi- 
tion, A can be tuned by an external magnetic %ld. The 
expressions for A and B can be found elsewhere.u The re- 
pulsive case is an idealization of the situation considered 
in Ref. ^ see Ref. [E§ 

We seek the solution in the form 



tp{x,y,z) 



n i 



;ikz 



c(x,y,z) 



+ b(x,y, z) e lkz , 



z < 
z > 



(2) 



where the amplitude a of the incident wave is fixed (real), 

j-2 i 2 

and the electron energy E — %r^. We assume that there 

w 2m 

is no current inflow along the screen (the only inflow into 
the system is from z = — oo). Thus only those solutions 
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satisfying the condition of zero inflow at 2=0, x, y — > ±oo 
will be considered. 

It is convenient to write (Jfy and (|l|) in dimensionless 
form by means of the definitions: x — \/2kx,y = \/2ky, 
z = kz, b = b/a, c = c/a. Insertion of (||) into Equation 
(0) for z ^ yields 



A ± & + i 9 52 & - % d s b = 0, 5 < 
A±c+ \d~ zi c + i dfc = 0, 2 > 0. 



(3) 



By using the continuity of the wave-function ip, one gets 
at z=Q: 



d s c - d s b + 2i (c - 1) = 2(q + /3|c| 2 ) c, 



(4) 



and A_l&=Aj_c. Here a=mA/(h 2 k) and (3=mBa 2 / (h 2 k). 
Eqs. (|^) and (0^ have spatially uniform solutions c=£ + i£ 
such that ( = —at; — /3£ 2 , |c| 2 = £, and 



/3 2 e + 2a/?r + (" + l)f - 1 = 0. 



(5) 



A straightforward analysis of this equation demonstrates 
that there is only one real root for a 2 < 3 and there are 
three real roots under the conditions: a 2 > 3, a[3 < 0, 
and P~ < (3 < (3+ with ^ = £ [t (a 2 - 3) 3 / 2 - a 3 - 9a] . 
Thus multiple solutions are expected for two cases: the 
barrier (a > 0) with attractive nonlinearity ((3 < 0) (case 
A); the quantum well (a < 0) with repulsive nonlinear- 
ity ((3 > 0) (case TV). Taking (3 as a control param- 
eter these solutions are depicted in Fig. [j] for different 
a. Notice that we obtain up to three coexisting uni- 
form solutions for different values of a: Z-shaped curves 
£(/?) (if a > \/3) and S'-shaped < a < 2) or loop- 
shaped (if a > 2) curves CO^)- At a=2 there is a cusp 
of the maximum of the £(/?) curve. The peaks in Fig. 
□(a) correspond to maxima of the transmission for which 
]c| 2 =£=l and (3— —a. Since (3 oc a 2 , multiple solutions 
exist on a certain interval of incident wave amplitudes 
for any strength of the nonlinearity B. The threshold 
values a=mA/kh 2 =±\/3 for multiplicity of uniform so- 
lutions can be achieved by varying the barrier height 
(well depth) and/or the energy of the incident wave. 
Three uniform solutions coalesce at the tricritical param- 
eter values: a =±\/3, (3 =T&V3/9, £o=3/4, ( 0=t ^/s/4. 
Hereafter we use the upper sign for case A and the lower 
sign for case 1Z. 

We shall perform now a small-amplitude perturbation 
analysis of Eqs. (||) and (Q) near the tricritical point. 
As a result we will find simple amplitude equations that 
will be solved in two particular cases of interest: (a) y- 
independent solutions, and (b) axisymmetric solutions. 

Let a = 0£q±5, (3 — /3oT7 with S > 0, 7 > and S, 7 <C 
1. We look for small nonuniform solutions: c = £ + i(, 
£, = £o+€i(x,y), ( = Co +(l(x,y), where £1 < f , Ci < 
£0 • The richest distinguished limit corresponds to having 
7 = Is + 0(^/ 2 ), Ci = O(VJ), 5, y = O^-Va), 
and z = 0(i5 _1 ). Inserting this Ansatz into Eqs. (^), the 
terms <9^& and d^c are 0(£ 5 / 2 ) and can be ignored when 
compared with the others, which are 0(<5 3 / 2 ). Inserting 
the result into (0), we find 
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FIG. 1. Real (a) and imaginary (b) parts of the transmit- 
ted amplitude c as functions of f3 for a uniform solution and 
different values of a. Only case A is shown. The case 1Z can 
be obtained by replacing a — > —a, /3 — » — /3. 



<9ij£l + dyy^l 



V3 



32 f 3 
27 ?1 

(f 7 -5) + 0(,5 5 / 2 ), 
Cr =±^/V3 + 0(<5), 



(6a) 
(6b) 



Notice that our Ansatz corresponds to weakly nonlinear 
perturbations of uniform solutions varying on a large spa- 
tial scale x=\[2kx — 0(S^ 1 ^ 2 ) ^> 1. The typical trans- 
verse length over which our solutions vary is thus much 
larger than the wavelength 1/fc. 

With the substitutions: £1 = §3 1/4 6 1/2 u, x = 
3 i/4^-i/2 X; y =3 i/4 (5 -i/2y ; Eq ^ Q can be written in 

the simpler form: dxxu + dyyu = u — 2u 3 + /1 + 0(<5), 
/i = 3- 1 / 4 <5- 3 / 2 (| 7 -(5) = O(l). We report here only the 
results (for y-independent solutions and for the axisym- 
metric case) corresponding to the most symmetrical situ- 
ation 7 = 1 8 (/i=0), where explicit formulas can be easily 
obtained. The results for the general non-symmetric case 
will be published elsewhere. 

a. Two- dimensional solutions depending on one 
transversal coordinate. If u=u(X) (two-dimensional so- 
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lutions of the full problem depending on only one 
transversal coordinate), the parameter- free equation 
dxxu = u — 2it 3 can be integrated once yielding the 



result (dxu) = u 



C. This equation admits 



nonuniform solutions satisfying the condition of zero flux 
as X — > ±oo only if C—0. In this case we obtain the so- 
lutions u = rj scch(AT — Xq), with r/=l for the soliton and 
r/=— 1 for the antisoliton. Next we choose Xo=0. Finally, 
using relation ( |6b[ ) our solution for the transmitted and 
reflected amplitudes on the screen will be 

c(x) = ^ eT" /e [1 + VSrjX sech(Ax) e*"/ 3 ] , (7a) 

b(x) = -~ e ±i7r/3 [1 - 3 V3?7A scch( Xx) e T ™ /6 } , (7b) 

where A=3~ 1//4 <5 1 / 2 and the upper and lower signs re- 
fer to the cases A and TZ, respectively. On the screen 
z=0 the difference between the cases A and TZ lies only 
in the phase of the wave function and not in the in- 
tensities: |c(i)| 2 = |[1 + V3rj\ sech(Ax)], \b{x)\ 2 = 
j[l — 3y3r]X sech(Ax)] (the small terms ~ A 2 are dropped 
consistently with our scaling). The different phase fac- 
tors give rise to drastic differences in the wave function 
outside the screen, as it will be shown below. We have 
also checked that the solutions are linearly stable when 
time evolution is considered subject to the boundary con- 
ditions discussed earlier. 

The amplitudes of the transmitted and reflected waves 
outside the screen can be found from (||) using as the 
boundary conditions their values at z=0 and ignoring 
the small terms due and <9;;6: 



c(x, z) 



1 



V47i 



c(x',0) e i ( £ - £ ') 2 /(42) d j/ j ( 8 ) 



and the expression for b(x, z) is the same once z is re- 
placed by —z. 

In our 2D problem the intensities of the reflected and 
transmitted waves are nonuniform in space in contrast 
to the ID problem, where they are constant. Denot- 
ing |5(i,z)| 2 = *?{l + * t (f,2)}, |6(i,i)| 2 - *°{1 + 
i& r (x,z)}, and using for c(i,0), b(x, 0) the soliton-type 
solutions (0), we obtain *?=3/4, *°=l/4. The nonuni- 
form parts of the intensities are given by 



cosf ^ rvjP + <f>f tt ] sech(Xx')dx' (9) 




A\z\ 



where a r =— 3, crt—y/3, and the terms ~ A 2 were dropped. 
The phase <f> in the argument of cosine is different for the 
attractive and repulsive nonlinearities: (case A): 4>t ~ 
j^tt, <j}+ — ~ ^7r; (case TZ): 4>t= ~J2 7T ' <K — ~T1 1T - 

Spatial distributions of the wave intensities are ob- 
tained by numerical integration of Eq. ([)]) and presented 
in Fig. H in terms of the scaled coordinates X = Xx, 
Z = AX z. The off-screen wave intensities are shown 
starting from certain nonzero values of Z. The wave 
intensity on the screen is shown as a thin strip in the 



middle of Figs. ||(a),(b). The results can be interpreted 
as follows. 
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FIG. 2. Density plots for the wave function intensities 
created by scattering off the self-induced nonuniform pattern 
on the screen at Z=0 with attractive (a) or repulsive (b) non- 
linearities. White (black) color corresponds to the maximum 
(minimum) of the intensity for the soliton solution on the 
screen (77 = 1), and vice versa for the antisoliton solution 
fa=-l). 

The uniform incident flow spontaneously produces a 
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nonuniform soliton-type pattern on the screen and then 
is diffracted by it due to the nonlinear feedback in the 
equations. In particular, for the soliton solution {r)=l) 
we observe local self-brightening of the transmitted wave 
with simultaneous local suppression of the reflected wave 
(Fig. [2]). The diffraction pattern is crucially determined 
by the value of the phase factor cf> ± . For the case A 
the transmitted wave is focused into a "beam" of higher 
intensity with a maximum outside the screen at Z w 1.7 
[Fig. [1(a)], whereas for the case 71 it is defocused and it 
"splits" into two "beams" [Fig. ||(b)]. Additional support 
for the importance of the phase factor is provided by the 
asymptotic behavior of the integral (^|) in the remote zone 
Z»l, 2X«Z 



~2 

! ^cos(^-- + 0%). (10) 
\z\ 4z 



4|z| 

t > for case 



In the transverse direction the local maxima (minima) 
of the intensities are determined by the condition 

<firt — 7rm > Tn = 0, ±1, . . .. For instance 
A, and the cosine in (|Io| ) reaches its maximum at x=0 
providing a transmitted "beam" along the axes. <fi7 < 
for case 1Z, and the cosine is largest on the parabola 
z = j^x 2 yielding two transmitted "beams" as in Fig. 
||(b) . The behavior of the reflected wave is also nontrivial: 
for case A the reflected pattern contains "split traces": 
the suppressed reflection forms the parabola z = -^x 2 
[Fig. ^(a)] , whereas for case TZ the reflection is suppressed 
within a single trace [Fig. §(b)]. It should be noted, how- 
ever, that for the antisoliton solution (r}=— 1) the maxima 
and minima are interchanged (with respect to the soli- 
ton solution), the "beams" become the suppressed traces 
and vice versa. Which type of solution (self-brightening 
or self-darkening of the transmission) will be realized in 
practice depends on additional conditions: type of the 
imperfections pinning the soliton, boundary conditions, 
past history, and so on. 

b. Two-dimensional axisymmetric solutions. For this 
geometry, u=u(R), R 2 = X 2 + Y 2 , and the parameter- 
free equation becomes Orru + -^9r u = U — 2u 3 . There 
are infinitely many solutions which satisfy our boundary 
conditions it(oo) = and ur(0)=ur(oc)=0. We denote 
them by S£ for the soliton case [u(0) > 0] and S~ for 
the antisoliton case [u(0) < 0], where the subscript n = 
0, 1, . . . is the number of zeros of as shown in Fig. |^. 
Then S~=—S+. 

In conclusion, spontaneous formation of spatial trans- 
verse patterns, which are quantum-mechanically coher- 
ent is expected to occur in semiconductor heterostruc- 
tures with a thin nonlinear layer. Self-diffraction of the 
electron wave on the transverse patterns gives rise to in- 
teresting phenomena like self-brightening or darkening of 
the transmitted wave, beam splitting, etc. 
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FIG. 3. Axisymmetric solutions S„ for the transmitted 
amplitude on the screen when n=l,2,3. Recall that Sn=—Sn- 
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